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Abstract: Motivated by the possibility of a string landscape, we reexamine tunneling of 
a scalar field across single/multiple barriers. Recent investigations have suggested modifi- 
cations to the usual picture of false vacuum decay that lead to efficient and rapid tunneling 
in the landscape when certain conditions are met. This can be due to stringy effects (e.g. 
tunneling via the DBI action), or by effects arising due to the presence of multiple vacua 
(e.g. resonance tunneling). In this paper we discuss both DBI tunneling and resonance 
tunneling. We provide a QFT treatment of resonance tunneling using the Schrodinger 
functional approach. We also show how DBI tunneling for supercritical barriers can natu- 
rally lead to conditions suitable for resonance tunneling. We argue using basic ideas from 
percolation theory that tunneling can be rapid in a landscape where a typical vacuum has 
multiple decay channels and discuss various cosmological implications. This rapidity vac- 
uum decay can happen even if there are no resonance/DBI tunneling enhancements, solely 
due to the presence of a large number of decay channels. Finally, we consider various ways 
of circumventing a recent no-go theorem for resonance tunneling in quantum field theory. 
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1. Introduction 

Various theoretical considerations have led to a picture of string theory with a gigantic 
number of solutions, referred to as the Landscape of vacua [1]. This picture was motivated 
in part by the construction of cosmological vacua in [2,3], and the dense "discretuum" in [4]. 
Although a detailed knowledge of the structure of the landscape is lacking, it is usually 
thought of as a complicated potential energy function in a multidimensional moduli space 
with a multitude of local minima separated by barriers of various heights. One of these 
minima corresponds to the vacuum where the universe exists today and confers it with 
its various characteristics (the cosmological constant, broken supersymmetry, the standard 
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model etc). Simple models of the landscape such as those discussed in [4,5] have different 
values of the cosmological constant and dimensionful couplings (such as Higgs mass) in 
the different vacua. The existence of a large number of vacua may shed some light on the 
cosmological constant problem [4]. 

Given an initial de Sitter vacuum with a string scale (or Planck scale) cosmological 
constant, successive tunneling events to vacua with lower vacuum energies relax the cos- 
mological constant. While this relaxation mechanism is happening, the universe at large 
scales eternally inflates and the global structure becomes very complicated [6], with vari- 
ous regions of spacetime inflating in various vacua. These regions are separated by domain 
walls. Eventually, some region of spacetime tunnels down to a metastable vacuum with a 
low cosmological constant, long lifetime and standard model fields. On the largest scales 
(far beyond scales which could ever be observed) the universe is not in our vacuum, and 
(assuming met ast ability) our vacuum will eventually catastrophically recede outside of our 
horizon, leaving a new and certainly hostile vacuum for us to observe. Tunneling is an 
important ingredient in this picture of the early universe. The universe samples the differ- 
ent sites of the landscape by some form of tunneling (two common well studied forms of 
tunneling being Coleman De Luccia and Hawking Moss instantons). 

Recent attempts at finding a cosmological measure for the landscape [7] employ tun- 
neling as the mechanism that allows the universe to sample various vacua. In fact, the 
occurrence of eternal inflation in the landscape itself depends on the mechanism of tunnel- 
ing.^ Tunneling mechanisms that have been widely studied in field theories [8-10] display 
an exponential suppression of the tunneling rate with an increase in the height of the barrier 
between the vacua. The expectation in the string landscape has been that the tunneling 
rates will be far suppressed compared to the rate of expansion of the universe during the 
de Sitter phase due to the generically string scale barriers. As a result, most regions of the 
universe eternally inflate, i.e. most regions are stuck in their various vacua with various 
positive vacuum energies. Certain regions of spacetime occasionally tunnel to regions of 
the landscape with a fiat potential conducive to slow roll inflation. Other regions that 
tunnel to vacua with negative cosmological constant are unstable to gravitational collapse 
(big crunch singularities) [9]. 

Because tunneling plays such a central role in the string landscape, it is important 
to fully understand it in such a setting. We assume from the start that transitions in 
the landscape have a four-dimensional effective field theory description^. In this paper 
we reexamine and expand the scope of tunneling in the context of the string landscape. 
In particular, we discuss two modes of tunneling that differ from the traditional view of 
tunneling in that they can lead to greatly enhanced transition rates even for tall barriers. 
These two modes are resonance tunneling [14] and DBI tunneling [15]. A crucial feature 
of these two modes of tunneling is that both of them make distant vacua accessible for 
tunneling. That is to say, in addition to the usual tunneling between neighboring vacua 
(those vacua which share a common barrier), tunneling can now occur at a relatively 

'^It is not clear if slow roll eternal inflation is generic in the string inspired inflationary scenarios. In 

some scenarios it is absent [11,12]. 

^ We thank Shanta de Alwis for emphasizing this point to us [13] . 
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unsuppressed rate between vacua that are separated by other vacua (i.e. vacua which are 
not immediate neighbors). We use the term Efficient Tunneling for this enhancement of 
tunnehng rate that can make transitions to far away vacua relatively unsuppressed [16]. 
This efficiency in tunneling (access to far away vacua) increases the number of decay 
channels from a given vacuum. When the number of decay channels becomes large enough, 
tunneling can become rapid (i.e. transitions from one vacuum to another can occur with 
order one tunneling rate). Efficient tunneling is a necessary (but not sufficient) condition 
for rapid tunneling. As we shall discuss in this paper, the string landscape can naturally 
provide an arena where efficient tunneling can also become rapid. 

This paper consists of the following parts. In Sec. |2| we discuss resonance tunneling 
using semiclassical methods. In the semiclassical approximation one can describe resonance 
tunneling in QFT models of the landscape. To this end, we use the methods developed by 
Ref [17-19] to study quantum tunneling in QFT. We discuss the importance of coherence 
for resonance tunneling to take place in a certain part of the landscape. Sec. ^ of this 
paper involves a discussion of DBI tunneling based on the analysis done in [15]. This 
provides a description of tunneling phenomena in the open string landscape where the 
dynamics is described by the DBI action. We shall see that both resonance tunneling and 
DBI tunneling involve major modifications to the Coleman-De Luccia result and lead to an 
enhancement of tunneling rates. Further, we shall see an interesting connection between 
resonance tunneling and DBI tunneling: DBI tunneling for supercritical barriers naturally 
leads to a condition necessary for resonance tunneling. In Sec. Q of this paper we use 
results from percolation theory to show how these efficient tunneling mechanisms, even if 
rare occurrences in the landscape, can still make a significant qualitative change to the 
picture. It is expected that a typical vacuum in the landscape will have a large number of 
decay channels. In flux compactifications, for example, a large number of decay channels 
is a result of many independent fluxes (cycles). The possibility of efficient tunneling, even 
if rare, together with the existence of a large number of flux directions can have nontrivial 
implications. In particular, in the presence of a large number of decay channels efficient 
tunneling can lead to rapid tunneling. Even though resonance/DBI tunneling offer novel 



ways of vacuum decay, they are not entirely necessary for rapid tunneling. In Sec. 4.1 



we show how the presence of a large number of decay channels, as expected in the flux 
landscape, can make even the usual vacuum decay a rapid process. 

In Sec. 1^ we discuss the various cosmological implications of our investigations. Here 
we summarize them: 



• The standard Coleman - De Luccia tunneling rate is responsible for the conventional 
view of the string landscape as an eternally inflating ensemble of meta-stable vacua. 
In this framework, tunneling to nearby vacua is exponentially suppressed, and tran- 
sitions to vacua with larger cosmological constant, or to vacua more distant in field 
space incur additional exponential suppression. Thus each site in the landscape has 
some degree of physical viability. While disconcerting for most purposes, the enor- 
mous size of the accessible landscape can be appealing for reasoning based on the 
anthropic principle or similar statistical methods. This picture may undergo striking 
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modification when efficient tunneling and the large dimensionality of the landscape 
are fully taken into account. In this new picture, two effects will alter any statistical 
framework. Firstly, certain decay channels are found to be vastly preferred over ev- 
ery other channel (sometimes over even those much closer in field space). Secondly, 
many vacua formerly thought to be meta-stable can become essentially unstable due 
to rapid tunneling. These two features can effectively remove large portions of the 
landscape, and may limit the domain of eternal inflation. 

• Fast tunneling needs to be shut off for vacuum sites with low cosmological constant 
(comparable to today's value). This could be achieved in two ways: by removing 
the possibility of efficient tunneling, or by reducing the number of nearest neighbors 
(dimensionality) in the field space. If this shutting-ofF of rapid tunneling is sudden, 
then this can be a natural explanation for the cosmological constant problem. The 
universe begins at some vacuum site with a large (string scale or Planck scale) vacuum 
energy. However, rapid tunneling makes the universe quickly decay down to a (non- 
eternal) inflationary vacuum site that ends with reheating and a low vacuum energy. 
In [16], RG group arguments are used to give an explanation for the necessary sudden 
shut-off of rapid tunneling. It is known that gravitational effects on vacuum decay 
can elegantly prevent [9] certain transitions to negative cosmological constant vacua, 
but they are of no help as a shut-off mechanism for the small domain wall tensions 
which coincide with rapid tunneling. Instead, one must assume another mechanism 
is able to induce the shut off of rapid tunneling which our current vacuum enjoys. 

• Fast tunneling gives a natural solution to the empty universe problem inherent in 
Brown-Tictclboim and Bousso-Polchinski scenarios. The empty universe problem 
arises due to the exponentially large lifetime of vacua from which the universe decays 
down to today's vacuum. The exponentially large time required for the decay dilutes 
away all the radiation/matter content of the universe. Bousso-Polchinski used the 
Hubble fluctuations to perch a classical infiaton field for subsequent slow roll inflation 
as a plausible solution to this problem. However, rapid tunneling provides a similar 
natural solution due to the short lifetime of vacua with high vacuum energies. This 
is because rapid tunneling can resemble classical inflation. 

• In some scenarios, rapid tunneling can make parts of the landscape accessible during 
the last 60 e-folds of inflation. This can lead to the formation of domain walls, the 
phenomenology of which has been discussed in [23] . 

2. Resonance Tunneling 

There are several examples where the quantum mechanical behavior of scalar fields plays 
an essential role during periods of cosmological inflation. Arguably the most familiar 
and important example are the inflaton's quantum fluctuations which lead to the density 

perturbations in our universe. Guth and Pi did a proper quantum mechanical treatment 
of the scalar field in [20] . The quantum treatment of a scalar field in the stringy landscape 
is further in that spirit. 
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In many situations, a quantum field theoretic (QFT) problem reduces to a simple 
quantum mechanical (QM) problem. For example, to understand aspects of the behavior 
of the infiaton field, it suffices under many circumstances to look at the analogous quantum 
mechanical problem [20] . The main complication that arises within QFT is the presence of 
an infinite number of Fourier modes. While the dynamics of QM is contained within the 
Schrodinger equation for the single particle wave function, the QFT problem is described 
by the Schrodinger functional equation, and the field is described by a wavefunctional. For 
a free field theory the wavefunctional factorizes to an infinite number of wavefunctions - 
one for each Fourier mode, and the QFT is then neatly described by an infinite number of 
Schrodinger equations. However, when the field theory is an interacting one (as it must be 
for tunneling to exist), these infinitely mimcroiis of Schrodinger equations become coupled 
with each other. Describing tunneling in a QFT might therefore appear to be a very 
complicated problem. Banks, Bender and Wu [17] extended the simple WKB tunneling 
picture in QM involving one degree of freedom to multiple degrees of freedom. Others 
[18, 19] extended the analysis to an infinite dimensional configuration space (QFT). The 
main idea underlying this extension from QM to QFT is that within the purview of the 
semiclassical approximation, tunneling under a barrier occurs in a tubular region within 
the infinite dimensional configuration space which can be described as a solution to a 
Euclidean classical equation of motion. This trajectory can be elegantly described using 
Euclidean path integrals, as in [8]. This enormous simplification of the problem under the 
semiclassical scheme makes the problem of tunneling tractable in QFT. In the following 
subsections we shall use the Schrodinger formalism to apply this basic idea to resonance 
tunneling within QFT. 

2.1 Resonance tunneling in quantum mechanics 

We first summarize resonance tunneling in quantum mechanics. Consider a particle moving 
in one dimension x under the influence of a potential V{x). Fig.(l) shows a potential with 
minima A, B and C separated by barriers. The wavefunction of the particle, V'(^)) can be 
obtained by solving the Schrodinger equation HiIj = Eip. In the classically allowed region, 
via the WKB approximation, this gives 



^l,r{x) ^ i— exp ^zfci j dxk{x)^ , (2.1) 



^fk(x) 

where k{x) = ^J'^{E — V{x)). The ± signs correspond to the left moving and the right 
moving solutions. A general solution in the classically allowed region will be a superposition 
il^{x) = aL'tpLix) + aRipR{x). Similarly, in the classically forbidden region the solutions are 

V'±(x)~— i=expf± / dxKix)] , (2.2) 

y/K{x) \ J J 



where k{x) = yJ^{V{x) - E), and the ± signs correspond to the growing and the decaying 
solutions. A general solution in the forbidden region will again be a linear superposition 
tp{x) = a+tp^{x) + a-^-{x). When a situation involves the propagation of a particle 
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in both the classically allowed and classically forbidden regions, the full solution has the 



form Eq.(2.1) in the allowed region and Eq. (|2.2D in the forbidden region. Since these both 
represent the same solution in different regions, the coefficients ol, or, and a_ are 
related by the connection formulas. For concreteness, let us consider the problem of the 
tunneling of the particle from vacuum A to vacuum B in Fig.(l). Let the coefficients of 
the components 'tpL,R in vacuum A be aL,R and the coefficients in vacuum B be Pl,r- 
These coefficients are related to each other by the following matrix expression which can 
be derived from the connection formulae (see, e.g., [21]). 





i(e-G-i 

9 + 9-1 




(2.3) 



where Q is given by 



~ 2 exp Q ^ ' dxy/2m{V{x) - E)^ , (2.4) 



xi and X2 being the classical turning points. To study tunneling of a wave packet from 
vacuum A to vacuum B, we keep only a right moving component in vacuum B (i.e. there 
is no particle incident on the barrier from the right), I3l = 0. The tunneling probability is 
given by 

TA^B = \^/ = 4(e + l,y\±, (2.5) 



9 / 92 



as is usually an exponentially suppressed quantity. 

Now we consider the possibility of resonance tunneling. We are interested in the 
tunneling of the particle from vacuum A to C via B. The two intervening barriers (between 
A and B and between B and C) contain the classically forbidden regions while B is classically 
allowed. One can again derive the relation between the coefficients ol^r (in vacuum A) 
and jL,R (in vacuum C) and express it in the following matrix form. 

if e + e-i i(e-e-i)\ 
" 2 1 -i{@-e-^) e + e-i j ^ 




-i($-$-i) $ + ^>-i 




(2.6) 



where 



and 



~ 2 exp Q y dx^/2m(yix) - E)^ (2.7) 
W= dx^/2m{E - V{x)) , (2.8) 

J xo. 
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with X3 and the turning points on the barrier between B and C. The transmission 
probabihty from A to C, via B, is now given by 

TA-.c = ^(^(e^ + ^^ + sin'^^ ■ (2-9) 

In general this will be a small quantity. However, when the so called resonance condition 
is satisfied, T^— >c can be enhanced. This happens when both W satisfies the "bound state 
condition" 

W = (^n+^^TT, n = 0, 1, 2, 3, ... (2.10) 

such that cos W = 0, and — > <I>, i.e. Ta-*b Tb^c- Under this resonance condition the 
transmission probability from vacuum A to C is greatly enhanced, Ta^c ~^ 1- Physically 
this enhancement is due to the constructive interference that occurs in the classically 
allowed region in vacuum B when the resonance condition (Eq.( |2.10| )) is satisfied. 

One can estimate the probability that the resonance condition is actually satisfied. 
We define AE to be the energy window where the resonance condition is roughly satisfied. 
This can be the range of E where Ta^c varies between 1/2 and 1 (i.e. it is the width at 
half maximum). Let Eq be the approximate energy difference between a typical pair of 
neighboring resonances. Then the probability for the occurrence of resonance is AE/Eq. 
Typically this is an exponentially small quantity, of the order of Ta^b and Tb-^c- More 
precisely [14], AE/Eq ~ {Ta^b +TB^c)/'2^■ 
2.2 Resonance tunneling in quantum field theory 

Tunneling in QFT has been well studied in the semiclassical regime using both the Eu- 
clidean path integral approach [8] and the functional Schrodinger equation approach [18,19]. 
The tunneling between two vacua is often well described by a WKB calculation. The tun- 
neling rate is generically suppressed as the height or width of the barrier between the two 
vacua is increased, and a typical tunneling rate is exponentially small. The tunneling pro- 
ceeds via bubble nucleation, and if the true vacuum interior is large enough, the bubble 
wall grows and quickly attains relativistic velocities. The picture might be more compli- 
cated when there are more features to the shape of the potential. Such considerations are 
especially important now in view of the string landscape. As recently discussed in [14,23], 
the presence of multiple vacua in the string landscape can give rise to various interesting 
effects arising due to efficient tunneling. One expects a large multiple vacua separated by 
high barriers to be a generic feature of the landscape (see [29]). In [14], a simple landscape 
is considered consisting of vacua A, B and C, with B being the intermediate between vacua 
A and C in the field space, as shown in Fig.(|^). When the resonance condition is satisfied, 
one expects the transmission probability to be of order unity to go from vacuum A to 
vacuum C via vacuum B. This means that any other vacua will have negligible branching 
fractions. The passage of the wave packet through vacuum B is quantum mechanical and 
leads to a constructive interference that is responsible for making the transmission proba- 
bility close to one. This, of course, can have nontrivial effects for the landscape and early 
cosmology (see e.g. [14,23]). 
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B 



C 



Figure 1: Vacua A, B and C separated by barriers in tiie landscape. 

In this section we begin a field theoretic analysis of resonance tunneling. The analysis 
in [14] is quantum mechanical and is based on solving the Schrodinger equation for a single 
particle. However, to understand the problem in QFT we must look at the functional 
Schrodinger equation and analyze the problem using WKB methods. The semiclassical 
approximation to the Euclidean path integral cannot capture the effect of resonance tun- 
neling. The Euclidean method captures the under-the-barrier effect between two vacua. 
However, for resonance tunneling, the motion in the intermediate classically allowed region 
is important. Although the Euclidean instanton technique is quite elegant in studying the 
problem of barrier penetration and false vacuum decay, the same problem can be studied 
using WKB methods in the functional Schrodinger equation formalism [18,19,22]. These 
references, in turn, build on the ideas developed by Banks, Bender and Wu in [17] to study 
quantum tunneling in dynamical systems with a large number of degrees of freedom. We 
shall adapt the methods spelt out in [18, 19] to study resonance tunneling in QFT. This 
method captures both the classically forbidden region as well as the classically allowed 
region. It is therefore especially suited to study the phenomenon of resonance tunneling 
which involves both forbidden and allowed regions in the field space in an essential way. 

The basic idea is based on constructing the WKB wavefunctions for the system in the 
multidimensional configuration space (corresponding to the various Fourier modes). This 
is simple for a system with a single degree of freedom (as in quantum mechanics for a 
single particle). The work done in [17 19] generalizes the WKB construction to systems 
with many degrees of freedom by arguing that under semiclassical approximation, barrier 
penetration happens mostly in small tubes within the configuration space around certain 
classical solutions, so that the WKB approximation is still essentially one dimensional. 
These classical paths correspond to classical solutions with Wick rotated (imaginary) time. 
The order h corrections can be calculated by considering small fluctuations about this 
classical trajectory. Ref [18] calculates WKB wavefunctions to the first two orders in h 
systematically. 

One begins with the full Hamiltonian H{(j)(x),(j)(x)) = H{(/){x), j6 / 5(l){x)), and writes 
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down the Schrodinger functional equation for the wave functional ^{(f){x)) of the field. 
Using the WKB approximation one may solve for this wave functional. The lowest order 
WKB approximation gives the equation of motion that determines the most probable es- 
cape path (MPEP) along which tunneling proceeds. The next order WKB gives corrections 
due to gaussian fluctuations around this path. In the remainder of this section we construct 
the WKB solutions in both the classically allowed and classically forbidden regions. 

2.3 Configuration space 

The configuration space of the field (l){x) is of infinite dimension. (In our 3 + 1 dimensions, 
X corresponds to the spatial directions.) We shall use A to parameterize a one-dimensional 
path in this configuration space. According to [17], in the semiclassical approximation, 
tunneling happens along a "most probable escape path" (MPEP) that can be parameterized 
by the variable A. This makes the problem tractable. The nearby paths will contribute 
quantum corrections that can be systematically calculated. The parameter A has nothing 
to do with the imaginary time that one uses in Euclidean techniques, it is just a parameter 
that describes the MPEP. 

A point in this configuration space is labeled by a parameter A. The coordinates 
of this point are given by the complete function {(/>(x; A)} = {(p{xi; A), 0(x2; A), ...}. Here 
{xi, X2, ...} span all of space and, hence, there is an continuous infinity of these. A trajectory 
is given by, for example, < A < 1, and as A continuously varies, so does the field 
configuration given by ^(x. A). A line element dr along the path (f){x, A) is given by 

{drf = J d\{dm? = {dXf J d\ ' . (2.11) 

So if the length along a path, r, is used to parameterize the path, then A = r, and 
J d^x {d<p/dr)'^ = 1. One can define the tangent vector along the path length dr and a 
vector orthogonal to the path ^(x, r) as 

dr ' 

= 6(1)- adr^ , (2.12) 
or 

where 5(p is an arbitrary variation, and a is chosen such that 

J d^xSd)^ (f^) = ' (2-13) 



kan = dr- 



which implies 

adr 



= (2.14, 



In general, there are an infinity of such orthogonal vectors at every point (^(x, r) of the 
trajectory. For every such direction n{y) we may define a unit vector as 
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where (t>n(y) {x-, f) depends on r, and the variable y denotes the coordinates of the direction 
n{y). The orthogonal variation can now be formally written as 



(2.16) 



Therefore, for a point (p{x) near the path defined by (p{x,r), there exists a value of r such 
that 



(p{x) = (p{x, r) + A(p±{x, r) 
where A(p_i_{x,r) is given by Eq.( |2.16| ). 

2.4 Functional Schrodinger equation 

The Hamiltonian for the scalar field (j){x) is given by 



.3 /li2 1 



(2.17) 



(2.18) 



where we assume that V{cj}) has the shape shown in Fig.(|T|). To quantize the field theory, 
one notes that the commutator of (p and (/) should be given by [(l){x), 4>{x')] = ihS^{x — x'), 
and so 4>{x) is replaced by the operator —ih5/5(j){x). The time independent functional 
Schrodinger equation is given by 



where the Hamiltonian operator is 



H 



(2.19) 



(2.20) 



Here ^{(f){x)) has the usual probabilistic interpretation, giving a measure of the probability 
of the occurrence of a field configuration (j){x). 

2.5 WKB approximation 

We assume a solution of the form ^'(0) = Ae'ii ^^'^\ where A is a constant. The Schrodinger 
equation then looks like 



/ 



d-^x - — 



h^^-¥[-ir) ^2^^^^ ^^^^^ 



e^^W = Ee^^W (2.21) 



Expanding S{(j)) in powers of h, S{<j)) = S'(o)('/') + ^'S'(i)(i;^) + and comparing terms with 
similar powers of h in the above equation leads to the following equations 



d-^x 



1 /SS(n)((t>)Y 1 



(2.22) 
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. ^'>g(o)(0) ^ 6S(o0)5S^^ 



(2.23) 



etc. 



Eqns.(2.22, 2 .231 ) represent an infinite set of coupled nonlinear equations in the infinite 
dimensional configuration space of the field (pix). To first approximation, however, it is 
possible to reduce these equations to a one dimensional equation. To this end, first define 
the potential ?7(0) related to V{(f)) as follows 



U{<t>) 



d'^x 



(2.24) 



The classically forbidden region is defined as the region where U ((/>) > E, and the classically 
allowed region is where U{(j)) < E. Let us first look at the WKB solutions for the classically 
forbidden region (region II of Fig.(^)). The analysis for the classically allowed region will 
be similar. The crux of the argument in [17] is that there is a trajectory in the configuration 
space of the field - the most probable escape path - perpendicular to which any variations 
of 5(0) vanish, and along which the variations of S'(o) are nonvanishing. If we use the 
parameter A to describe this one dimensional path in the infinite dimensional configuration 
space, then the most probable escape path is given by (poipc, A), and 5S(j^-)/5(f) is nonvanishing 
only along this path. If such a path exists, then we can convert the hopelessly complicated 
set of Eqns(|]2|, |2^...) to a very simple one dimensional equation. This can be done as 
follows. The most probable escape path satisfies 



5S, 



(0) 



5S, 



(0) 



\(l>o{x,\) 







(2.25) 



Further, C(A) is related to dS^Q-j/dX, and can be eliminated. Note that along the most 
probable escape path one can write 



(0) 



d^x- 



5S, 



(0) 



ix,X) 



dX 



(2.26) 



so 



d^x 



dX 



(2.27) 



Combining Eq (|2.27| ) and Eq (|2.25| ) with Eq( p.22D , one gets the WKB equation 




dX 



dX 



U{(t>o{x,X))-E . 



(2.28) 
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Finally, using the definition of the path length dr in the field configuration space as given 
by Eq( 2.11| ), one reduces the field theory WKB approximation to the one dimensional 

1 /dS, ^ 2 



^ - U{r) - E, (2.29) 



2 \ dr 

where U{r) is the potential along the most probable escape path given by 

U{r) = j d'x Q(V0o(x,r))2 + y(</.o(x,r))') . (2.30) 



Similarly, Eq. (|2.23| ) for 5(1) becomes the one dimensional equation 

_d'^S(n\ _ dS(n\ dS, 



Q^2 Qj. Qj. 

The familiar solutions to Eq(2.29) and Eq.(2.31) are 



(0) ,-(0)-(y_^^ 



5(0) = i / dr'^2[U{T')-E], 



5(1) = -\n^2[U{r)-E]. (2.32) 

The WKB wavefunctional effectively reduces to a WKB wavefunction and can now be 
written as 



^nm^, r)) = exp (^-i dr'y/2[U{r') - E] + , (2.33) 



A_ 

{2[U{r~E]y 



where the II subscript corresponds to the classically forbidden region II in Fig.(||). One 
can also determine the equation for the most probable escape path. This path corresponds 
to the Hamilton-Jacobi variational equation 

pr 

(55(0) =^ dr'^y2[U{r')- E] = (2.34) 
Jo 

for arbitrary variations (j)± away from (j)Q. The resulting path is then given by the Euler- 
Lagrange equation 

+ V</.o(x,A) — = 0, (2.35) 

where A is related to r by the relation 

- = ^2[U{X)-E] . (2.36) 

Note that (a) A is a real parameter, as we are concerned with the tunneling region where 
U > E, and (b) A is simply a parameter we introduced to solve the equation of motion 
and, in particular has no relation with the Euclidean time. If a solution to Eq.( p.35 ) exists, 
then it can be used as the classical trajectory that describes the most probable escape path 
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along which tunnehng occurs. An infinite dimensional problem thus gets reduced to one 
described by a single parameter A (or, equally well, by r). 

One can similarly construct the WKB solutions in the classically allowed region of the 
field configuration space. The WKB wavefunction for the classically allowed region is given 
by 

where the subscript "c" in (j)c{x,s) means that this classical solution corresponds to a 
classically allowed trajectory. Here s parameterizes the path. The classical orbit can be 
determined by solving the Euler-Lagrange equation 

5!^_vV=(..,) + ^^^^^ = 0. (2.38) 

where rj is defined as 

^ = ^2[E-U{v)] . (2.39) 

2.6 Corrections to WKB solutions due to neighboring paths 

Although the main contribution to the WKB approximation comes from the classical path 
(either Euclidean or Lorentzian), there will be additional contributions to it from neigh- 
boring paths. In the case of tunneling, as discussed originally in [17], the tunneling takes 
place in a tubular region that encloses the classical Euclidean path and also contains vari- 
ous neighboring paths. These neighboring paths give a calculable correction to the WKB 
solutions. We shall simply quote the results here. The interested reader should refer to [19] 
for the details. 

To describe the corrections, one must define coordinates that are orthogonal to the 
tangent along the classical path in the field configuration space. These orthogonal coordi- 
nates n{y) will be used to describe the nearby paths. Given a classical path (pQ{x,r), one 
can describe a nearby path (/)(x, r) as (from Eq.(^l6|)) 



(j){x, r) = (j)o{x, r) + A(j)±{x, r) , (2.40) 
where the perpendicular displacement (with respect to the (/)o(x,r)) is given by 

A0±(x,r) = j d^y4)^(^y){x,r)n{y) , (2.41) 

with (l)n{y){x, r) being the unit vector in the orthogonal direction n{y). Here y describes the 
coordinates along the orthogonal direction. Thus infinitesimal variations of an arbitrary 
path (j){x, r) at constant value of r are given by 

x)\r = (i)niy){x,r)Sn{y) . (2.42) 
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Further, let us define the second variations of the action in the classically forbidden and 
allowed regions as respectively 



K{y,y';r) = 
K{y,y';s) 



d"^xid^X2 



'^'5(0) (0) 



{Xi)5(j)±{x2) 



(0)1 



6(t)±{xi)6(l)±{x2) 



\M^,r) X(i>niy){xi,r)(i)n{y'){x2,r), 
U,(x,s) X^n{y)ixi, S)^n{y'){x2, s) ■ (2.43) 



Including the effect of neighboring paths, the wavefunctional for the classically forbidden 
region around the classical tunneling path is given by 



^{^{x,r)) 



A ( 1 
-77 exp — - 

{2[U{r)-E]f'^ V ^ 



dr'yj2[U{r')- E] + 



- d y d y'K{y,y',r)6n{y)6n{y') + .. 



(2.44) 



Note that the paths away from the classical path give exponentially suppressed contribu- 
tions, and this justifies the leading order expectation that tunneling takes place in a tubular 
region around the classical solution. Similarly, the wavefunctional for the classically allowed 
region around the classical orbit is given by 



^i(p{x,s)) 



A 

777 exp , ^ 



ds'^/2[E -U{s')] + 



- d% d'y'K{y,y',s)6n{y)6n{y') + 



(2.45) 



m(A) 



d^x 



The above WKB solutions can be written in a slightly different form which serve better 
for studying the matching conditions c.f. Sec. (2. 7). One defines a "mass" parameter as 

Note that the mass parameter is always positive, m(A) > 0. Note that from Eq.( 2.1lj ), 
dr/dX = \J m{X). Now one can define the WKB solutions both in the tunneling region and 
the allowed region via the same parameter A as 



A 



i2[UiX) - E]f" 
1 
2 



exp 



in the forbidden region and 



^{4>{x,x)) 



A 



{2[E-U{X)]fl^ 
1 
2 



d)^ yj2m{X')[U{X') - E] + 
d^y I d^y'K{y,y',r)6n{y)6n{y') + . 



dX'y/2m{X')[E -U{X')] + 



(2.47) 



expl- 



d^y / d^y'K{y,y',s)6n{y)5n{y') + ... 



(2.48) 



in the allowed region, respectively. 
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Figure 2: The potential U{X) as a function of the parameter A. 



2.7 Matching conditions 

The greatest simphfication that has happened in the above analysis is that tunneling in 
the infinite dimensional configuration space has been reduced, to the leading order, to a 
quantum mechanics problem in one dimension. The effect of nearby paths appears first at 
next to leading order. Let us first consider the problem of tunneling across one barrier, and 
then the case of tunneling across two barriers (and the condition for resonance tunneling): 

• Tunneling across a single barrier: As in Fig.(||), consider three regions, regions I, 
II and III. The energy E is such that regions I and III are classically allowed whereas 
region II is classically forbidden. As in [14], let an and q_l be the coefficients of the 
right -moving and the left-moving components, respectively, in region I. Similarly, let 
and be the coefficients in region III. The two turning points are given by Ai 
and A2. Then the WKB connection formulae imply 





(2.49) 



where Q is given by 



e 



2exp(- 



dX 



(2.50) 



and where Ueff{X) is the result of integrating out the fluctuations (denoted in the 
previous section by 5n) around the tunneling path (i.e., the Euclidean classical solu- 
tion). Ueff{X) = U{X) + F{X), where F{X) contains the one-loop effect that appears 
when one integrates out the fluctuations denoted by 5n. To evaluate F we use 



exp 



F{T')dT' 



det' 



'1/2 



(2.51) 
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U(X) 




Figure 3: A potential with multiple barriers and minima. Resonance tunneling can happen as a 
result of constructive interference. 



where the prime on det' implies, as usual, that the zero mode is to be excluded while 
evaluating the determinant. The variable r is related to A by ^ = \J 2[jT(Xy^ '' ^'^'^ 
we write ti = r(Ai) etc. The difference between the simple WKB tunneling formula 
in quantum mechanics and the one in quantum field theory is the inclusion of the 
gradient term |(V(/>)^ and the one- loop effect F. The transmission probability from 
region I to region III is then given by 



Ti- 



>iii 



4(G + 



1^ 
9' 



-2 



4 

02 



(2.52) 



Resonance Tunneling: We consider the tunneling across two barriers, Fig(^). Re- 
gions I, III, and V are classically allowed, whereas regions II and IV are classically 
forbidden. Again, one can follow through the analysis of [14], keeping in mind that 
one should also include the gradient term and the one-loop effect. Let and 7^ be 
the right-moving and the left-moving coefficients in the region V. 



Then we have 




where 



and 



$ ~ 2 exp 



i(e-G-i 

+ 0-1 





A4 
A3 



dX'^J2miX')[UeffiX')-E] 



W= dX'J2m{X')[E-Ueff{X')] . 



(2.53) 



(2.54) 



(2.55) 
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The transmission probability from region I to region V via region III is now given by 




r,^v = 4 ( ( 9$ + ^ ) cos^ + + sin2 . (2.56) 

Again, although this resembles the quantum mechanical result in Sec.(p.l|), the crucial 
difference in field theory is that one must consider the effective potential Ueff{4>) in 
all the WKB results and not the potential Vifp)- Now one can easily see the resonance 
condition 

W=(^n + ^^7T, (2.57) 

where n = 0,1, 2, 3, .... 



3. Branes and Tunneling 

Here we discuss the role branes play in vacuum transitions, necessitated by their central 
role in realistic backgrounds in flux compactifications. The primary ingredient in stable 
string compactifications are wrapped branes, and the charges they impart on the various 
cycles of the 

internal geometry. A subset of these branes are mobile in the sense that they are 
point-like in the six-dimensional bulk. Both the point-like and wrapped branes contribute 
a significant fraction of the landscape [24, 25] , although their dynamics is quite different 
than the closed string moduli which control the topological and geometric data of the 
vacua^. The brane dynamics is decribed by the DEI action which has a nonstandard kinetic 
term. Decay processes involving the brane, therefore, are described by DBI tunneling 
which, as we shall see, can differ significantly from the Coleman-De Luccia tunneling [15]. 
Before describing DBI tunneling, we will briefly explain other vacuum transition processes 
involving branes. 

The Brown- Teitelboim instanton [26] describes a tunneling event whereby a volume- 
form flux is changed discretely. Charge conservation dictates that the volume form be 
discontinuous across the domain wall by an amount equal to the brane charge, /i. The 
domain wall in this case is the charged brane itself, and its tension is r. Because the brane 
is of zero thickness, the decay rate through such a process is exactly described by the 
thin wall limit of Coleman-De Luccia. The volume-form flux which couples to the brane 
contributes an amount ^ 

V=-*{F4A*F4) (3.1) 

to the vacuum energy, and so the energy difference between the two vacua is given by 

e = ^,y^-^^l\ (3.2) 
where V- is the contribution to the true vacuum energy due to F^. 



^In F-theory, however, the closed and the open string landscape can be treated in a unified manner. 



A more exotic form of branc-induccd vacuum decay is the Kachru-Pearson-Verlinde 
instanton [27]. In this scenario, a flux stabilized, "doubly contractible" spherical brane 
tunnels by slipping around a larger sphere, and the potential barrier height is given by 
the brane tension times its spatial volume. As expected, a higher potential causes greater 
suppression of the decay rate. 

3.1 DBI tunneling 

DBI tunneling is distinct from the previous two scenarios in that no quantum numbers 
need be changed. Instead, a single mobile brane simply tunnels between meta-stable loci 
within the internal geometry. What is unusual about this instanton, is that the decay rate 
increases when the potential barrier height is elevated to near the (five dimensional) BT 
threshold. This effect is due to the decreasing bubble wall tension, to which the decay rate 
is exponentially sensitive. To illustrate this, we consider a mobile D3 (or anti-D3) brane in 
a warped bulk geometery. Here and below we follow [15] closely. The action for the brane 
in a warped background is given by ^ 



where 1/ f{(f)) is the warped brane tension. The potential F(^) has two important sources: 

(1) any five-form flux component which is incompatible with the supcrsymmetry preserved 
by the brane, and (2) threshold corrections from Kahler stabilizing D7 branes. 

We would like to compute the rate per unit volume for the three-brane to transition 
between two local minima of the potential V((p). The rate for this process is given by the 
path integral over the Euclidean action 



We are using the semi-classical approximation, where the field corresponding to the classi- 
cal Euclidean solution (pd dominates the path integral. Proper normalization is achieved by 
subtracting the Euclidean action of the constant field in the false vacuum, and we assume 
the fluctuation determinant sets the overall coefficient to order one in fundamental units. 
We call the field (t)ci "the bounce," and it is the unique non-constant solution to the equa- 
tions of motion with lowest Euclidean action which asymptotically approaches the false 
vacuum at Euclidean infinity. Such a solution has been shown [30] (under reasonable as- 
sumptions) to be a maximally symmetric interpolating solution, i.e. its location breaks four 
generators, but preserves the remaining 0(4) symmetry of the underlying (Wick rotated) 
Poincare group. This allows us to write the Euclidean action as 




(3.3) 




(3.4) 



—Se[4'ci]-Se14>+] 



(3.5) 




(3.6) 



''The essential features of DBI tunneling do not depend on the warp factor. 
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where dots represent derivatives with respect to p. 

We focus on two neighboring minima of l^((/>), e.g. A and B of Fig. (1). A vacuum 
transition is much hke a first order phase transition, and begins with the formation of a 
bubble of true vacuum within the false vacuum. The bubbles will tunnel quantum mechan- 
ically to a critical radius and then expand classically at essentially the speed of light, 
converting the vacuum energy difference into kinetic energy of the bubble wall. The high 
degree of symmetry has reduced our computation to the one-dimensional motion of a parti- 
cle in a potential — y ((/>), albeit with an unusual kinetic term. In the language of Sec.(^), 
the solution to this one dimensional equation of motion gives the most probably escape 
path (MPEP). In our analysis, we will use conjugate momentum, ti^, and Hamiltonian, H, 
which are defined with respect to L^^'^°'^"^'^y(27r^/j^). 



L = ^,^ = ^ + y , (3.7) 



2^ canonical ^ 

dL ■ , ^ 

TT^= —r- = (t)-f , (3.8) 



H = TT^<P-L = - V , (3.9) 



where 



As TT(f) — > f{4>)~^^'^ the field velocity ^ ±00. This means that when one plots the in- 
stanton in p-cj) coordinates, <j) is not a single valued function of p. This failure of coordinates 
is due to the fact that the world line in Euclidean space continues to curve in response to 
the forces imparted by the potentials, and p need not be always increasing along the world 
line. One encounters a similar situation with the instanton that describes Schwinger 
pair production in the presence of an external electromagnetic field where the Euclidean 
classical solution is double valued as a function of the Euclidean time. This, of course, is 
not surprising as the relativistic action for a charged particle in an external electromagnetic 
field is the one dimensional version of the DBI action that we are considering here. The 
Euclidean equation of motion is given by 

vr^ = -— -3-. (3.11) 

Notice that 

dH = --7r^# (3.12) 
P 

which determines the rate of non-conservation of H. 

Smoothness of the brane at p = requires that (j) begin (and end) with zero velocity 
(7 = 1), so dH integrates to the energy difference between the two vacua. A fraction of 
this energy is matter inside thick-walled bubbles. For thin-walled bubbles, the energy is 
purely vacuum energy, causing the well known "empty universe problem." 
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Thin walled bubbles occur whenever Si ^ e///, where e is the difference in vacuum 
energy, fi is the mass of (j) (in either vacuum), and Si is the tension of the domain wall. 
Then we can consider the energy nonconservation to be negligible (i.e. dH = upto 0{e) 
in Eg. ( p. 12] )). The mass fi sets the scale for the bubble wall thickness, and the ratio Si/e 
determines the bubble radius c.f. Eqn.( 3T^ ). In the inverted potential, cj) sits exponentially 



close to the true vacuum, </>_, for a long time (p ^ 1/^) before rolling off and continuing 
quickly to the top of the false vacuum, </)+. The thin wall bounce is approximately a 
step function. A first integral of the equation of motion Eq.( 3.1l| ) is given by writing 
H + 0{e) = E in Eq.(|!|), yielding 

7(</., vr^) = 1 - /(0) (F(0) + E + 0{e)) , (3.13) 

where E = — V+, the value of the potential at the false vacuum. This is solved for the 
momentum to give 

vr0 = yWHWWW)) , (3.14) 

where we have absorbed E -\-0{e) into the potential by setting both minima of Vq equal to 
zero. The only obstacle to modeling friction as a conservative force is that it is momentum 
dependent rather than position dependent. But if the momentum can be shown to be 
purely a function of position, as is done with Eq. (|3.14 ), the obstacle disappears, and we 



may write an effective (but trajectory-dependent) potential which completely describes the 
forces, and thus motion. Defining Vq is simply absorbing the friction into the potential, a 
feat possible because the bounce (/>(/)) is a monotonic map (in the sense that its inverse is 
a well defined function). The term of order e is linear in cf) to first approximation, and so 
we can write 

Vo{<t>) = V{4>) - V{<f>+) + (3.15) 
From this we can find the solution by inverting the function 

1 - f{4>)Vo{(t>) 



(3.16) 



A refinement of Vo(</>) can be made iteratively by solving for H{(j)) using this solution, and 
plugging it into Eq.( |3.13| ). 

These solutions determine the bubble nucleation rate through the difference in Eu- 
clidean action B = Se[(I)] — SE[(j)+]- Notice only the interior volume and bubble wall 
contribute to B. 



B = -27T^ r p'^edp + 27T^f V^o(2-/Fo)# , 
Jo J<j>- 

= -^e + 2^V5i , (3.17) 
where p is the radius of the bubble. We may integrate Eg. ( ^.121) to solve for 

p = 3Si/e. (3.18) 
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The tunneling rate per unit volume is given by 

The bubble wall tension is 5*1 = J d(/)y^Vo(2 — /Vq), so the tunneling rate is large when 
/Vo — > 2 over much of the potential. This tunneling rate is always greater than the 
Coleman - De Luccia rate in the same potential. Although in general this does not imply 
that the tunneling is rapid, but if fVo is close to 2, the rate can be substantially greater 
than Coleman-De Luccia. 

3.2 The Wrinkle 

As one might expect from a Euclidean solution, there is a generic multi-valuedness of (j){p). 
The complete range of 7 is given by 

-1 < 7 < 1 =^ < f{(P)Vo{(P) < 2 . (3.20) 

A typical solution would look something like Fig.(^). It is clear why the solutions for 
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Figure 4: Two solutions for the instanton trajectory, (dashed) for a low barrier (/Vo < 1) and 
(solid) for a high barrier (1 < fVo < 2). 



large potentials is "wrinkled," i.e. folds over and becomes multivalued. This is because the 
extrinsic curvature of the solution is determined by the gradient of the potential. This is 
the Euclidean version of -F = ma. While finite acceleration prevents an initially time- like 
world line from wrinkling, no such effect occurs in Euclidean space (as is evident from the 
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S instanton that exists for the Schwinger pair production for charged particles in external 
electromagnetic field). 

One might worry that the presence of the wrinkle (and the consequent large value 
of the second derivative of (f)) implies a breakdown of the DBI description of the brane. 
However, the correct measure of how big the a' corrections are is the extrinsic curvature of 
the brane and not just the second derivative of cf) with respect to p. The DBI action can be 
trusted whenever curvatures are low [28]. In the thin wall regime, the extrinsic curvature 
of the solution ^(p) in the warped background is given by the covariant derivative of a 
unit one-form n normal to the curve ^(p), whose pullback is traced over using the induced 
metric. 



The bulk metric is given by the Euclidean warped throat: ds'^ = a' f~2dp'^ + a''^d(j)'^. For 
the validity of the DBI description K{(p) should be smaller than the string scale 1/y/a'. 
This limits how steep the potential can be while still trusting the DBI description. 

3.3 Including the effects of gravity 

Gravity complicates the calculation in a number of ways. The main effect is that the 
background is no longer fiat, although we still assume that an 0(4) symmetry is preserved 
in the maximally symmetric, uniformly curved background. Another effect is that the 
friction felt by the rolling (p is now dynamical, and coupled to the inflaton Hamiltonian. 
Although this coupling makes the general solution rather difficult to approximate, the thin 
wall limit remains explicitly calculable. Here we simply quote the results [15]. 

The decay rate in the presence of gravity (with Mp the reduced Planck mass) is given 



K{(f)) = {V^h^)h''^diX^djX' 




(3.21) 



by 



V/V ~ e 




(3.22) 



where 




v=v. 



(3.23) 



5wall = 27r2p35, . 



(3.24) 



Here 



35i Mp 



(3.25) 



P = 



^e^l + ISl {V+ + V.) + l,St/Ml 



and 




(3.26) 



as before 
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3.4 Efficiency 

The Brown- Teitelboim instanton is responsible for neutralizing any potential which exceeds 
a certain critical threshold analogous to the pair-creation threshold. One might expect this 
to lead to a large population of vacua which, despite having high potential barriers, are 
below the BT threshold. It is precisely such an ensemble of vacua in which DBI tunneling 
is much more rapid than the standard CDL result, but still slow enough to be within the 
semi-classical region of validity. One would expect that the truly rapid case (i.e. where 
tunneling is so rapid that the term meta-stability is no longer operative) exists in two 
distinct limits. The first is simply the small barrier behavior near a critical point, i.e. 
where a first order transition nears second order. The second is where resonant effects 
become large, causing the barrier to resemble a reflectionless potential (we shall see in the 
next section how the DBI behavior can lead to this). ^ 

3.5 Resonant DBI and the Moonscape 

Two interesting features appear when DBI tunneling is combined with super-critical po- 
tential barriers, i.e. those for which /Vq > 2 over any part of the barrier. The first is 
that the BT mechanism allows for nucleation of D3 bubbles of various topologies. For 
relatively low barriers, the only allowed topology consistent with our 0(4) invariant ansatz 
is X S"^, but for higher potentials an S'^ bubble will be dominant. (The requires a 
higher potential since it must occur near p = 0, where friction is large.) The second, and 
more striking feature is that mobile D3's will see the barrier as effectively two separate 
barriers of equal height (see Fig.(|5|)), since the interior supercritical region is classically 
allowed. This doubling of vacua naturally leads to a setting where resonance conditions 
can potentially be satisfied. It is this possibility which enables the transition to be much 
more dominant than the closed D3 bubble solutions, although Euclidean methods cannot 
show this. Additionally, this provides a mechanism for pairwise similar barriers, causing a 
statistical enhancement for the resonance condition to be satisfied. The resulting picture 
is that the effective maximum barrier height corresponds to Vq = 1//, since all higher 
potentials contribute less to the domain wall tension density. The reduced tension can be 
described by the CDL-effective potential 

Vo^'''' = Vb(l - fVo/2) , (3.27) 

where the CDL (canonical) kinetic term is assumed. In the following subsection we will 
show that the correct way to calculate the domain wall tension, and thus nucleation rate 
for supercritical barriers will be to use the above effective potential Vq"^^ with a standard 
kinetic term, via the methods of section |2|. Thus high 'mountains' in the open string 
landscape are turned into 'craters,' since the only suppression occurs in the window < 
/Vo < 2. The open string landscape is really a moonscape. Because the supercritical 

third, more speculative limit exists where the potential rises so steeply as to bypass the intermediate- 
height regime, i.e. a step function. It should be noted that this third limit only exists far outside the range 
of validity of the DBI action, and we would have to assume that the unsuppressed a' corrections do not 
conspire to prevent these rapid transitions. We will not consider this possibility here. 
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barrier regions contribute nothing to the domain wall tension, long distances in field space 
are just as accessible as short ones, a criteria which can lead to percolation. 



V 








Figure 5: The dashed curve represents the actual potential energy Vq for a supercritical barrier, 
and the solid curve is the effective (non-relativistic) potential Vif^^ seen by the field. The many 
craters appearing from high open string barriers make up "the moonscape." 



3.6 Craters and resonance tunneling 

As we have just shown, the DBI action leads to an enhanced tunneling rate which can be 
traced to the fact that the domain wall tension is corrected to 5i = / d(py/Vo{2 — /Vq)- 



Equivalently, the potential can be reduced to a CDL-effective potential given by Eq. (|3.27| ). 
In fact, this equivalence is not just limited to D-branes. In this section, we show how 
this equivalence can be easily seen in the case of a point particle obeying the Klein-Gordon 
equation. For moderately high barriers, the effective potential resembles a crater, and when 
the barrier height exceeds the BT threshold, the crater interior will be a classically allowed 
region, as in Fig.(|^. As pointed out in [15], it is in such cases where no interpolating 
Euclidean solutions exist, although we expect the transition rate to be non-zero. This has 
been referred to as 'Klein tunneling' [39], which we illustrate in Fig.(|6|). In the absence 
of resonant effects, the closed (BT) instanton is dominant, but satisfying the resonance 
condition will make the tunneling solution greatly favored. To calculate the transition 
rate, we will resort to a one-dimensional example of large barrier penetration by solving 
the Klein- Gordon wave equation. 

Consider a particle of unit charge and mass m interacting with an electromagnetic 
potential. For simplicity, we assume no warping. The action is given by 

S = -^J dxdt {'D^D''^ + m^"^^') , (3.28) 
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p-p 








Figure 6: Two solutions for the supercritical (/Vq > 2) case, an x bubble, and the "Klein- 
Tunneling" solution with a potentially resonant interior. The solutions differ only by initial condi- 
tions, since the transition must begin and end in the two vacua. The x D3 bubble dominates 
over the D3 transition, unless the resonance condition is satisfied. The S'^ D3 bubble occurs far 
below, near p = 0. 



where 




(3.29) 



The resulting Klein-Gordon equation D^D^^ = m?^! can be solved using separation 
of variables. Here we consider a single positive frequency mode of energy ui. The choice 
^ _ Q^'^tip(^x) removes the gauge symmetry At ^ At + dta, and so we shall replace V with 
the gauge fixed Vq. The time independent Klein-Gordon equation is then 

We have written this equation in a form which allows for manifest mapping of the Klein- 
Gordon equation to the Schrodinger equation. In so doing we have defined the Schrodinger 
effective potential V(^'^^ = Vq {1 — Vq / {2uj)) in complete analogy with Eqn.( ^^ ). It is 
now clear how we may use the non-relativistic formalism to describe relativistic resonance 
tunneling. In particular, DBI tunneling through a single barrier will be sensitive to resonant 
effects whenever the barrier is higher than twice the local D3 tension. 



4. Percolation in the Landscape: Efficient Tunneling and Preferred Paths 

The interesting consequence of both resonance and DBI tunneling is that when certain 
conditions are met, tunneling rates become highly enhanced. This effect can be traced 
to a reduced domain wall tension. In the case of resonance tunneling, tunneling becomes 
insensitive to the height and the width of the barrier (provided the resonance condition is 
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exactly satisfied), and tunneling proceeds "classically," i.e. without any suppression. The 
time scale of tunneling will be the natural scale of the problem. In cosmological scenarios, 
this time scale will the Hubble time. Similarly, in DBI tunneling if the right conditions 
are met (i.e. 1 < /Vq), tunneling becomes much more efficient. This is exponentially 
faster than the time scale afforded by usual tunneling, although still meta-stable in the 
controlled regime. When /Vq > 2, DBI tunneling naturally allows for resonant tunneling 



(Sec.( |3.4| ) and Sec.(p.5[)). In fact, as we later explain in Sec. 4T, even in the absence of 
any resonance / DBI tunneling, the enormous number of decay channels expected in a 
flux landscape can make vacuum decay a rapid phenomenon if certain conditions are met. 
Such efficient tunneling will have nontrivial consequences for cosmology, especially in the 
context of the string landscape which provides a huge number of vacuum sites and decay 
channels, and where the DBI action describes a large fraction of the transitions. In [14], 
Tye has suggested there might be a possible solution to the cosmological constant prob- 
lem using such efficient tunneling, [23] uses efficient tunneling to extract some interesting 
domain wall phenomenology from the last few e- folds of inflation (discussed in Sec.(^)). 
The crucial question, of course, is how frequently the conditions conducive to efficient tun- 
neling are satisfied in the landscape. If it is extremely rare to find conditions for efficient 
tunneling, then one might question the efficacy of efficient tunneling. In particular, mere 
isolated instances of efficient tunneling in the landscape will not alter the overall tunneling 
dynamics across the landscape. 



In this section we shall argue, along the lines of [14], that even if the relative probabil- 
ity of efficient tunneling is low, the string landscape might have the right characteristics to 
make use of it. Two facts will be crucial for this : (i) a large number of field (or flux) di- 
rections present in the string landscape that lead to a large number of decay channels, and 
(ii) the possibility to resonance tunnel between somewhat far away sites of the landscape 
(and, as discussed in Sec.(|3.4D and Sec. (3. 5), DBI tunneling can itself become resonant). 



The second condition is not entirely necessary. We explain in Sec. 4.1 how even usual 
vacuum decay can become rapid in the landscape provided certain conditions are satisfied. 
To quantify our discussion, we shall use some basic results from Percolation Theory [31] 
and explore the possibility of finding paths across the landscape that are continuously con- 
nected by channels of efficient tunneling. Such paths will make the landscape percolate. 



The basic idea underlying percolation theory, and its relevance to the landscape and 
efficient tunneling, can be illustrated by a simple example. Consider a two-dimensional 
grid of electrical resistors. The lattice points are the points where four resistors join (see 
Fig(|^)). Let this be a large enough grid so that the linear size of the grid is much larger 
than the linear size of a unit cell of the lattice. In other words, the number of resistors 
involved is a large number A^. Considering any two neighboring points on this lattice, let 
p be the probability of the resistance between the two points being zero, and so 1 — p is 
the probability that it is infinite. The question we ask is : starting from an initial value of 
p = (i.e., all resistors nonconducting), we slowly increase the value of p (i.e. randomly 
convert some of the infinite resistors to zero resistance so that it can perfectly conduct 
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Figure 7: The small circles represent lattice sites. Solid lines are paths of infinite resistance. 
Dashed segments offer zero resistance. For conduction, all the segments should have zero resistance, 
i.e. critical probability — 1. 



Figure 8: Bethe lattice with the number of connections to any lattice site given by z = 3. 

electricity). At what critical value of p, the critical probability does the grid become 
conducting? 

At this critical value, there will appear one continuous connected line of zero resistance 
resistors that will let the current pass through. The two dimensional grid will separate into 
two non conducting regions, separated by a conducting line. In other words, the grid 
will 'percolate. For two-dimensions, it is known that pc = 1/2 [31]. Of course, for one- 
dimension, pc = 1 (i.e., all the resistors have to be conducting; any nonconducting resistor 
will disrupt the flow of current (see Fig(|^). The value of pc for higher dimensional grids 



z = 3 
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Figure 9: A square lattice of resistors. All resistors have infinite value. No current flows. 



is not analytically known. However, when the number of dimensions is high, the answer 
may be easily approximated. For large number of dimensions, a hypercubic lattice in d- 
dimensions behaves like a Bethe lattice (Fig.(^)) with 2d nearest neighbors, for which the 
value of pc is exactly known. For a Bethe lattice with z nearest neighbors, the critical value 
is 

Pc = ^ (4.1) 

z — I 

Since a higher dimensional hypercubic lattice (in d-dimensions, for large values of d), can be 
approximated as Bethe lattice with 2d nearest neighbors, the value of the critical probability 
is 

Pc = — (4.2) 
^ 2d-l ^ ^ 

Note that as the number of dimensions increases, pc decreases. 

To apply the above idea to the string landscape and efficient tunneling, let us consider 
a hypercubic lattice of string vacua motivated by the Bousso-Polchinski flux landscape [4] ^ 
This approach lets us talk more generally without committing to any particular realization 
of the landscape. All that is required is individual vaccum sites connected by decay channels 
- a situation easily represented by a Bethe lattice. We shall be solely interested in efficient 
tunneling and shall consider the usual tunneling that takes an exponentially large time 
as "no tunneling". As we are interested in the efficacy of efficient tunneling, such an 
approximation is valid. Let p be the probability of the occurrence of efficient tunneling. 



^Although our picture is motivated by the Bousso-Polchinski landscape, the mechanism we describe is 
much more general without any reference to the particular realization of the landscape. In the Bethe lattice, 
seen in Fig.(^), the lattice points would represent the vacua and the lines connecting them would represent 
the decay channels. 
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Figure 10: Solid lines are infinite resistors. Dashed lines are zero resistance connections. The 
number of dashed lines are too few to make the grid conducting, i.e., p < Pc- 



Typically we expect p to be a small number. For resonance tunneling p will be exponentially 
suppressed (as discussed in Sec.( ^.lD ). For DBI tunneling p can be much higher as it 



naturally gives a resonance condition for supercritical barriers (Sees. (3^ and |3.5| )). 

The question we ask is : for a given value of p, what is the dimensionality d of the 
hypercubic lattice (i.e., the number of flux directions, or in a more general set up simply 
the number of tunneling channels) for which percolation happens. Naively, as we saw in 
our discussion of Bethe lattice, the answer is l/{2d — 1). Hence, for a typical value of d in 
the 100 to 1000 range (corresponding to a hundred to a thousand flux directions / decay 
channels), the critical probability is of the order ~ 0.01 — 0.001. This might make DBI 
tunneling percolate if p is in this range and if the DBI tunneling is highly enhanced. On 
the other hand, resonance tunneling has an exponentially suppressed value of p (Sec. (2.1)). 
However, as we will now argue the answer can be exponentially better for resonance tun- 
neling, enough to make it percolate. 



There is one important difference between the electrical conductivity of the hypercubic 
lattice of resistors and quantum tunneling in the landscape. In the grid of resistors, a 
current can only flow from a given lattice point to a neighboring lattice points. However, 
in the landscape quantum tunneling can occur between non-nearest-neighboring sites as 
well. The number of sites to which a wavepacket can potentially quantum tunnel to can be 
estimated. Let $g be the length in the field space over which quantum coherence can be 
maintained (i.e. the length of the field space over which the phases of the wavefunctional 
have definite relationships). For example, during infiation, for light fields it is natural to 
assume that H as they have a natural spread over a Hubble scale. Let $0 be the 

typical spacing between neighboring sites in the landscape. Then the number of sites to 
which a wavepacket can potentially resonance tunnel to, starting from some initial site, can 
be approximated as ~ {^g/^o)'^, where d is the number of directions about the particular 
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Figure 11: At critical probability pc, a continuous curve of zero resistance appears through the 
lattice. The lattice becomes conducting. Note the contnuous dashed line (zero resistance path) 
from lattice site A to B. In 2-dimensions, Pc = 1/2. 



landscape site. This can be an exponentially large number. For the choice of = 3<I>o 
(i.e. coherence is maintained over three lattice sites) and d = 100, we get = 3^^^. Now, 
while using Eq. (|4.1| ) to find the critical probability for percolation to happen, one should 
use z ~ iV. 

Then we have 

Pc-^^- (tA ' (4-3) 



N \<^>_ 

For the coherence length spanning just a few landscape lattice spacings ^q, a very small 
value of Pc can be achieved provided d is large. 

The scenario discussed in [14] assumes an exponentially small value of p, the probability 
of the occurrence of efficient tunneling. The probability of hitting resonance, p, is shown to 



be of the order of the usual (nonresonance) tunneling probability (as discussed in Sec. (2.1)). 
That is, p is exponentially suppressed. For example, a typical of tunneling probability 
(using the thin wall Coleman-De Luccia result) in the Bousso-Polchinski scenario is ~ e"^''" 
[4]. Assuming such a small value of p, we see from Eq.( |4.3| ) that percolation can still happen 
as the critical probabilty p^ can itself be low. There can be a continuous path connecting 
string landscape sites across which resonance tunneling can happen. The existence of 
such continuous paths allowing high transition rates between vacua can have important 
cosmological consequence. Most interestingly, there can be no eternal inflation from such 
vacua due to the large number of decay channels, and the decay rate not being so highly 
suppressed. 
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Figure 12: The dots are the vacua. Starting from vacuum A, quantum tunneling can potentially 
happen to all vacua enclosed in the dotted coherence sphere. Solid directed lines imply that non- 
resonant tunneling can happen along the arrow. The dashed directed line implies that resonance 
tunneling can happen from A to C via B. is the lattice spacing. 

4.1 Rapid tunneling via the Chebyshev metric, and toy landscapes 

There is another means by which the large dimensionality of the landscape, in combination 
with efficient tunneling, can yield qualitatively different behavior. When a large number 
of decay channels are present, the total decay rate is written 

r ~ ^ iV, exp {-27T,''stj2e^) , (4.4) 

where Nq is the number of channels of a given rate (or tension, if we fix e). For definiteness, 
we will assume that the tension Siq is of the same order as the vacuum energy difference e. 
Each decay channel is therefore very slow, of order exp(— 277r^/2) fa 10~^^. Qualitatively 
new behavior occurs when the number of accessible final states is large enough to cancel 
this factor. Let us assume that the field configuration space has 300 dimensions, and is a 
regular square lattice. We will imagine that the domain wall tension is determined through 
some measure of distance in the field configuration space, i.e. the single domain walls 
associated with neighboring vacua are of lower tension than the composite domain walls 
which bridge to more distant vacua. Pictorially, Nq is the surface area of a shell in the field 
configuration space which surrounds the initial vacuum. 

By reformulating the domain wall tension as a function of distance in field space, the 
physics is contained in the choice of metric on the (discrete) space of vacua. We will 
consider several different metrics, all being discrete versions of norms. These metrics 
represent distinct toy landscapes, where the binding energy of various domain walls (which 
determines the tension of the composite domain wall) is of the simple form of an I'p norm. 
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Figure 13: To apply percolation theory, tunneling from vacuum A to tlie other vacua inside 
the coherence sphere in the hypercubic lattice of the previous figure is better represented using a 
Bethe lattice. The number of other lattice points connected to a given lattice site is z, which is 
approximately the number of vacua inside the coherence sphere. 

The three cases we consider will be the unbound (or marginally bound) case p = 1, the 
BPS case p = 2, and p = oo, known as the "Chebyshev" case. 

If one assumes that all of the composite domain walls are marginally bound, i.e. have 
no binding energy, then 

Si [qi,q2, -qsoo] = Si x {\qi\ + \q2\ + ... + |g30o|) , (4.5) 

and the shell will be a sphere defined by the "taxicab metric" on the field configuration 
space. If on the other hand the domain walls bind together {d la BPS) such that 

5*1 [qi,q2, ■■■qsoo] = Si x ^ ql + ql + ... + ql^^ , (4.6) 

then the shells will be spheres defined using the standard metric of R'^'^''. We will think 
of the BPS formula for the domain wall tension Eq. (|4.(]| ) as a realization of moderately 
efficient tunneling. 

The more interesting "Chebyshev" example would have a tension 

Si [qi,q2, ...gaoo] = X maxdg'il, |g'2|, kaool) , (4.7) 

in which case the shells are hyper-cubes. What makes this version of efficient tunneling 
interesting is that the volume remains large even when the number of dimensions is large. 
In a field configuration space of such large dimension as the landscape (say 300), there is a 



-32- 



dramatic enhancement for the highly efficient Chebyshev case. This is because the surface 
area of a "taxicab" sphere of radius p is quite small in spaces of large dimensionality. 



Fo/(5^^^) = 2(2g)^«7300! (4.8) 

taxi 



The standard sphere enjoys the larger surface area 



yoZ(5'^^) = 3007r^^^g'^7l50! (4.9) 

Euclid 

but a truly large volume sphere is the Chebyshev sphere: 

FoZ(S299) = 300 X 2300g299 (4 ^0) 

Chebyshev 

A simple two-dimensional potential which exhibits Chebyshev domain walls is given by 

V{cj)i,cj>2) = sm\cj)i) + sm\<j>2) - 2 sm\cj)i) sm\cj)2) . (4.11) 

With such a metric Nq {2q)^'^^, and it is easy to overcome the large suppression in 
the Coleman - De Luccia tunneling rate. If we assume that the individual domain wall 
satisfies Sf/e^ = 1 then the total decay rate is given by 

few 

T ^^Nq exp (-277r2g4/2) (4.12) 

q=l 

^ 10^^ (4.13) 

This rate is so large that the system would easily percolate. (Here we use the term 
'percolate' to mean the nucleated bubbles rapidly coalesce, rather than in the sense of 
percolation theory used previously.) Thus instead of eternal inflation in such regions of the 
landscape, the outcome would more likely resemble a thermalized domain wall foam. While 
this toy model is highly speculative, it shows how the large dimensionality of the landscape 
can conspire with efficient tunneling to transform eternal inflation into percolation. 



5. Discussion 

In this paper we have discussed two novel tunneling mechanisms in the string landscape - 
resonance tunneling and DBI tunneling. Both lead to an enhancement of the tunneling rate 
provided certain conditions are met. We also saw the interesting connection between the 
two: DBI tunneling for a supercritical barrier can naturally become resonant. Resonance 
tunneling makes tunneling to far way vacua in the landscape relatively unsuppressed, i.e. 
far away vacua becomes accessible and tunneling becomes efficient. Further, the existence 
of a large number of decay channels in the landscape can convert this efficient tunneling 
into a rapid one. Efficient tunneling is, therefore, a prerequisite for rapid tunneling. Such 
rapid modes of tunneling, where the tunneling rates are not exponentially slow compared 
to the other time scales in the problem, can have nontrivial consequences for the evolution 
of various regions of the universe populating the vacua in the landscape. As discussed 
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in Sec-d^ using basic ideas from percolation theory, even if the occurrence of efficient 
tunnehng is rare in the landscape, the existence of a large number of decay channels can 
make efficient tunneling rapid. 

In the Bousso-Polchinski scenario, the number of non-collapsing decay channels quickly 
decreases as one approaches the surface of the zero cosmological sphere (the Bousso- 
Polchinski sphere) in the lattice of flux vacua. As pointed out in [14,16], This reduction 
in the number of decay channels for vacua with very low vacuum energy (one of which 
corresponds to our current vacuum with the low value of the dark energy) can explain why 
inspite of efficient tunneling, the universe is stable today. Paramount to the success of 
this mechanism is a superselection effect which forbids transition to negative cosmological 
constant vacuua. The reduction in the number of decay channels makes efficient tunneling 
an improbable occurrence. In other words, away from the Bousso-Polchinski sphere the 
large number of decay channels can make efficient tunneling natural, whereas close to the 
Bousso-Polchinski sphere efficient tunneling can get shut off due to a drop in the number of 
decay channels. If this transition is sudden, then this picture gives a natural explanation for 
a small cosmological constant. It is important that the shutting-off of rapid tunneling close 
to the Bousso-Polchinski sphere be a sudden one, otherwise the intermediate region (where 
tunneling is slow, i.e. where metastability exists) can support eternal inflation. [16] uses 
RG group arguments - drawing a parallel with the well studied metal-insulator transition 
in condensed matter - to show that such a shut-off of rapid tunneling (which is analogous 
to the metal- insulator phase transition) is indeed very sharp. The argument rests on the 
properties of the beta function and the RG group flow in the vicinity of the region where 
the rapid tunneling occurs. One could give similar arguments using percolation theory 
drawing on the connections between percolation and renormalization group (see [31]). A 
superselection sector boundary at zero cosmological constant has been hypothesized by 
Banks [32]. 

If efficient and rapid tunneling is a generic feature of the landscape due to a large 
number of decay channels (away from the Bousso-Polchinski sphere), then it can interfere 
with the occurence of eternal inflation. Given any vacuum with a large number of decay 
channels, rapid tunneling will make it decay on a time scale that is exponentially smaller 
than that required for the occurence of eternal inflation. As efficient and rapid tunnelings 
constrain eternal inflation, they also constrain the general utility of anthropic reasoning. 
Anthropic reasoning assumes that the vacua in the landscape are generically populated 
due to eternal inflation. However, rapid tunneling makes eternal inflation impossible and 
can make a large number of landscape vacua unpopulated. This could dramatically alter 
the anthropic reasoning that is often adopted in landscape studies. 

As we saw in Sec.(|3|), DBI tunneling makes tunneling enhanced when the barriers 
are of stringy heights. This is due to the Schwinger effect that becomes operative due to 
the square root kinetic term (similar effect is seen of charged relativistic particles in strong 
electromagnetic field where the action has a square root kinetic term) . After taking account 
of this relativistic effect, the tunneling rate is the same as the Coleman-De Luccia tunneling 
but with a new potential Vif^^{cl}) which is related to the original potential Vo{(j)) by the 
relation Vq"^^ = Vo(l — /Vo/2), f{cl)) being the warp factor. What is also interesting is that 
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for supercritical barriers (/Vq > 2), DBI tunneling naturally leads to a resonance condition. 
Such a connection was certainly not obvious at the outset. Clearly, this connection between 
DBI tunneling and resonance tunneling deserves further investigation. In Sec.( |3.6| ), we start 
with the Klein-Gordon wave equation for a relativistic point particle (a zero-dimensional 
membrane) with a certain external potential V. We then show that the problem can 
be mapped to a Schrodinger equation with a modified potential V{1 — V/2). This is an 
alternate way of justifying the simple recipe Vq^^^ = Vb(l — /^o/2). 

Even though the details of resonance tunneling in quantum field theory are complicated 
compared to quantum mechanics because of the infinite dimensional configuration space, in 
effect all that happens is the constructive interference in the intermediate vacuum/ vacua. 
In order to meaningfully talk about definite phase relationships over distances in field space 
which are involved in the constructive interference, coherence has to be maintained in the 
field disturbances over that field distance. In Sec. (Q) we have represented this coherence 
length as $g and there is a corresponding notion of the coherent sphere. During inflation 
light scalar fields have coherence over a field distance ^ H. This can allow for the 
possibility of constructive interference over distances in the field space that are < H. 

An important issue in resonance tunneling is what happens to the conventional picture 
of vacuum decay via bubble nucleation. In [8], the vacuum decay is a first order phase 
transition that occurs via the nucleation of various bubbles of various sizes. The bubbles 
are of subcritical, critical and supercritical sizes. Subcritical sizes are too small and due 
to the excess surface energy in the bubble wall they collapse. Critical bubbles have an 
exact cancellation of the surface energy and the drop in the volume energy due to the 
bubble nucleation. Supercritical bubbles have a large enough size that it is favorable for 
the bubble to grow. This picture is bound to change for the decay of the false vacuum 
via resonance tunneling [14] . Resonance tunneling involves the quantum state of the field 
in the intermediate vacuum - the state corresponding to the constructive interference. [14] 
has argued that subcritical bubbles would nucleate with field values corresponding to the 
intermediate vacuum and these small bubbles would collide with each other leaving behind 
the final vacuum and radiation. The slow roll infiation phase could be immediately preceded 
by such a phase transition that released radiation. Possible observational signature of such 
a thermal phase preceding the slow roll inflation has been investigated in [34]. 

In a recent paper [43], the mechanism by which resonance tunneling might proceed in 
quantum fleld theory was investigated. The authors reported a no-go theorem that pre- 
cludes the occurence of resonance tunneling as a first order phase transition. The authors 
, however, have interpreted their result as a general no-go theorem for resonance tunnel- 
ing in quantum fleld theory. We have addressed this no-go theorem in Appendix The 
complete understanding of what should replace the usual bubble nucleation for resonance 
tunneling is lacking. However, as we explain in the appendix, the usual framework that 
involves an 0(4) symmetric bounce solution will breakdown. 

During resonance tunneling, the state of the fleld in the intermediate vacuum is quan- 
tum mechanical. The motion of the field in this region, in the WKB approximation, is given 
by a wavefunctional such as Eq.( p.48 ), an oscillatory solution in the field space. If the in- 
termediate vacuum resembles a quadratic well, then these classical solutions correspond to 
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the ground or the excited states of the field in the simple harmonic potential. Solutions 
that are very excited correspond to classical like oscillatory solutions of the field in this 
potential. The field starts at the initial false vacuum, undergoes classical like oscillations 
in the intermediate vacuum, and then exits at the final vacuum. It is possible that the 
oscillatory motion in the intermediate vacuum can lead to parametric resonance [35]. This 
will lead to particle production. The backreaction of the particles produced can destroy 
the resonance condition. It would be interesting to quantify the conditions under which 
the backreaction is large enough to interfere with the resonance tunneling. 

Efficient and rapid tunneling can have interesting phenomenological consequences [40] 
^. It will be interesting to do a detailed analysis such as [42] with rapid tunneling taken 
into account. One such scenario involving domain walls has been studied in [23]. Typically, 
since tunneling is assumed to be an exponentially slow process, any signs of nearby vacua 
are relegated to the far past when the universe was eternally infiating. This makes any 
possibilities of finding observational signatures of the string landscape impossible even 
in principle. However, efficient tunneling that becomes rapid can change this picture and 
provide some handle on nearby vacua. In particular, while the infiaton field sources inflation 
of the universe, there could be other light fields that have typical de Sitter spreads in their 
vacuum expectation values (a fiuctuation of the order of a Hubble mass scale during one 
Hubble time). This spread in the vevs can lead these light fields to explore near by vacua 
(the coherent sphere discussed in Sec.(^) is related to this spread in the vev of the light 
fields) and display rapid tunneling if the conditions are conducive. Assuming that rapid 
tunneling happens during the last 60 e-folds of infiation, different Hubble patches of the 
postinflationary universe could end up in different vacua. These vacua will be separated by 
domain walls whose tensions depend on the structure of the local landscape that is sampled 
by the efficient tunneling. Depending on the distribution of the heights of the barriers, 
there could be a distribution in the tensions of the domain walls resulting from efficient 
tunneling. As discussed in [23], requiring the domain walls to not pose any cosmological 
problems can lead to constraints on the structure of the local landscape (of course, this 
is under the assumption that efficient tunneling happens during the last 60 efolds). The 
vacua on either side of the domain wall will generically differ in their vacuum energies. If 
this split in vacuum energies is large enough, the domain walls will experience a pressure 
gradient and will attain relativistic speeds very quickly. This can help get rid of domain 
walls. Requiring the domain walls to go away before big bang nucleosynthesis (BBN) 
requires the split in vacuum energies to be greater than MeV^ and a lower bound on the 
domain wall energy scale of the order 10 TeV^. However, if the split in the vacuum energies 
is less than MeV^, the walls will stay around and eventually collapse into primordial black 
holes. The mass of the black hole will depend on the domain wall tension. A spread in 
the domain wall tensions will lead to a spread in the masses of the primordial black holes. 
Domain walls with higher tensions will lead to black holes with lower masses which, in 
turn, evaporate first. When the black hole evaporates it converts the energy which was 
initially in the domain wall into radiation. This leads to a reheating phase of the universe. 

''See also [41] and references therein. 
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In [23] it is shown that, starting with domain wahs of various tensions that are formed right 
after inflation, as soon a domain wall of a particular tension becomes the dominant form of 
energy in the universe, it collapses to black holes. These black holes eventually evaporate 
away. The lighter black holes will evaporate first. This can lead to multiple phases of 
reheating in the early universe. Requiring the reheat temperature to be high enough for 
successful BBN requires the wall tensions to be above 10~^Mp. Since the domain walls 
(and the black holes to which they collapse) become the dominant form of energy in the 
universe at the time of the collapse, it is not possible for these primordial black holes to 
be dark matter candidates. 

Finally, it is plausible that both modes of rapid tunneling might serve as alternate 
solutions to the empty universe problem, which one encounters in trying to explain the small 
cosmological constant using relaxation mechanisms [4,26]. The final stage of the sequence of 
sampling of different sites of the landscape should correspond to the very low cosmological 
constant that we observe today. However, tunneling typically takes an exponentially large 
time and, when implemented to solve the cosmological constant problem (a la Brown- 
Tietelboim, Bousso-Polchinski) has the associated problem of the Empty Universe. The 
original Brown-Tietelboim mechanism had two problems : the gap problem and the empty 
universe problem. The Bousso-Polchinski proposal succeeded in solving the gap problem 
due to the large number of directions that exist in the flux landscape, and gives plausible 
scenarios for solving the empty universe problem. The exponentially large time required 
for tunneling from one vacuum site to another dilutes away the matter and radiation 
components of the universe and leads to an empty universe - this is the empty universe 
problem. However, as we have discussed in this paper, tunneling need not always be an 
exponentially slow process. While tunneling usually takes an exponentially large time to 
take place, when certain conditions tunneling rates can be highly enhanced. Further, the 
large number of directions present in the landscape lead to a large number of decay channels, 
some of which might lead to rapid tunneling (as discussed in Sec.(^). In scenarios where 
such efficient and rapid tunneling might be generic, tunneling rate can be unsuppressed 
and the cold empty universe may be avoided. 

A relevant issue while discussing preferred paths (percolation) in the landscape is that 
of the initial condition ^. While in our analysis of resonance tunneling in Secs.(^?l], 2.2) we 



solve a time independent Schrodinger equation, the problem of tunneling is a dynamical one. 
One should prepare a wave packet of the particle/field and solve the corresponding time 
dependent Schrodinger equation. This is a standard exercise in quantum mechanics, one 
starts with an ansatz 'ip{x, t) = 0(x)e*'^* for the wave packet and reduces the time dependent 
equation to a time independent one. The results obtained from the time independent 
analysis are then recovered. Perfect resonance condition is satisfied only in certain bands 
AE of the energy of the wave packet. The wavepacket itself has a spread in energy 5E. For 
the complete passage of the entire wave packet through the resonance band the requirement 
is 6E < AE. This will ensure that even without a large number of decay channels, the 
efficient tunneling is rapid. What decides the width of the wave packet in the landscape? 



*We thank Puneet Batra for detailed discussions regarding the importance of the initial conditions. 
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Presumably this question is tied up with the question of initial conditions which has been 
much debated [37]. Simple applications of quantum cosmology lead to very different results 
about the initial condition of the universe [36] . As discussed in [38] this might be a problem 
of staying within the minisuperspace approximation. Going beyond the minisuperspace 
scheme can solve some pathologies of the wavefunction of the universe. Luckily, in our 
picture we do not need to know the width of the wavepacket. Resonance ensures that the 
tunneling is efficient. And the presence of a large number of decay channels can make the 
efficient tunneling rapid. 

In our analysis of resonance tunneling in QFT, we have not included gravity. In our 
semiclassical analysis we have assumed a fixed flat spacetime background. Admittedly 
gravity should be included while discussing tunneling in the string landscape. Inclusion 
of gravity will include a gravitational Hamiltonian along with the usual matter Hamilto- 
nian {H = Hm + Hg) in the Schrodinger equation. This Schrodingcr equation is simply 
the Wheeler-De Witt equation (H^ = 0) for the wavefunctional of the gravitational and 
the matter degrees of freedom expressing the diffeomorphism invariance. Now * is the 
wavefunction of the universe. However, we do not expect gravity to change the qualitative 
picture. For vacuum sites with a large vacuum energy, the effect of gravity is not impor- 
tant [9] . Gravity becomes important for vacuum sites with small vacuum energy by making 
them stable. 
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A. Appendix: The thin wall approximation and the landscape 

We have used the thin wall approximation throughout this paper and one might be con- 
cerned regarding the thin wall approximation when the vacua and their energy differences 
span large intervals ^. The energy difference between two neighboring vacua, e, will in 
general not be small in the landscape. However, the thin wall approximation can still be 
used when e is not exceedingly small. In this appendix we show that the thin wall formula 

^We thank Alexander Westphal for detailed discussions about the robustness of the thin wall approxi- 
mation. 
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Figure 14: The red (continuous) line shows the potential V{(j)) with a metastable vacuum at 0+ 
and a true vacuum at (f>-. The blue (dashed) line is a hypothetical inverted quadratic potential 
V = Vo — 5M^</'^ chosen so as to fit the actual potential V{(j)) as closely as possible in the field range 
between 0+ and (f>^. 



has a wide range of applicability by deriving a rough estimate of the upper bound on e 
(in units of the height of the barrier Harrier ) that still allows a thin wall description. In 
particular, we shall show that one can stay within the thin wall approximation and still 
have e ~ Vbarrier- We are interested in an order of magnitude estimate. Consider a poten- 
tial V{(p) with metastable minimum at 0+ and a true minimum at (p- (see the red curve 
Fig.(|l^). The energy difference between the two vacua is e = V{4>+) — y ((/>_). We assume 
that the potential V{(p) has a shape that allows a thin wall approximation for small e. We 
investigate what happens as e is increased starting with a small value well within the thin 
wall regime. 

For small e, the thin wall picture of bubble nucleation holds, i.e. the wall of the bubble 
is thin (compared to the radius of the bubble) and at the moment of the nucleation of the 
bubble, the field configuration outside the bubble is (p = (p-^. and inside the bubble it is 
(j) = (j)-. This configuration (the bounce) is obtained by solving the Euclidean equation of 
motion (using p for the radial Euclidean coordinate) 

^ + -?-m = (A.1) 
p dp 



which describes the motion of the field in an inverted potential —V{(j)). The boundary 

dp ' 



conditions for the bounce solution are : ^(/> = 0) = 0, and (t){oo) = (p+. In the thin wall 



approximation the friction term, is neglected. The bounce solution (j){p) begins near 
the true vacuum (j) = (p- at p = 0, stays very close to the true vacuum till it reaches the 
bubble wall at the radius p = R, then it quickly crosses the wall (the time over which it 
does so is given by the inverse curvature of the potential at its top) and goes very close 
to the false vacuum (p = (p^ and reaches this value at p = oo. A configuration such as 
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Figure 15: A thin wall bounce configuration. The radius R of the bubble wall is much greater 
than the thickness of the wall (the range of p over which the configuration changes over from 
to 



this is shown in Fig. (|T5|) . Over the wall region the field quickly rolls over the valley of the 
inverted potential —V{(j)) (i.e. the top of the potential V{(j))). The thickness of the wall is 
given by (5 ~ /x""*^ where /i^ = |V/op(i;^)| (i.e., /i is the curvature of the top of the potential 
barrier). The thin wall condition is therefore : R = 35i/e >> 



Now we start increasing e. There will be a gradual transition from the thin wall 
regime to the thick wall regime. For large values of e the bounce will take the thick wall 
configuration described in Fig. (16). At p = oo, the field will sit at its false vacuum value 

As we start moving towards the center at p = 0, the field will tend towards the true 
vacuum (/>_!_ but will never attain that value. Instead for all e above a certain value, the field 
configuration at p = will be (/)(0) = where can be uniquely determined by solving 
the Euclidean equation of motion (Eq.(A.l) with the boundary conditions (/)(p = oo) = 04-, 
and ^(p = 0) = At the nucleation of the (now thick walled) bubble, the field at 
the center of the bubble will be at = (/>*. Post nucleation, as the bubble grows, will 
classically roll from (j)^ towards (pj^. 

Next we find a rough estimate for (p^, , and the corresponding value of energy difference 
e* = y((/)+) — V{4)if). The thin wall approximation breaks down much before V{(j)-[.) — 
V{(j)-) = e^:. We shall determine this value etw shortly. The reader should note the 
difference between etw (above which thin wall is a bad approximation), and e* (above 
which the decay rate becomes very insensitive to To find e* we adopt the following 
strategy. We fit the potential F((/>) with an inverted quadratic potential V = Vq — ^p^(/>^ 
(blue line in Fig.(|l^)) For such a model potential the Euclidean equation of motion is 



^"intuitively it is clear that once e is very large, the decay rate will not depend on e any more. There 
is no way for the field sitting at the false vacuum 0+ to gauge the depth of the true vacuum on the other 
side of the barrier prior to tunneling. In other words, once e exceeds some critical value which is large, the 
decay rate will simply become insensitive to the exact location of the true vacuum. 

^^This analysis is therefore applicable to those potentials that can be fitted with an inverted quadratic 
potential in the region between the metastable and the true vacua. One could do a similar analysis for the 
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Figure 16: A thick wall bounce configuration. At p = oo, (p — (pj^. As p decreases, the field value 
tends towards 4>- but never quite attains that value. At p = (center of the bubble), <j) = <j)^. The 
value of 0* is obtained by solving the equation of motion for the bounce with appropriate boundary 
conditions, and satisfies the condition 0+ for our potential. The radius R of the bubble 

wall is the same as the thickness of the wall. 



linear and exactly solvable. Let be the point on this inverted quadratic potential that 
is the closest to the true vacuum of the actual potential V{(j)) that we are interested in. 
We assume that our original potential V{(j)) has a shape that allows the existence of such 



a very close to Now we solve the Euclidean equation of motion (Eq.(AJ)) using 
this inverted quadratic potential 

-i + --f+;uV = (A.2) 

where z is a Euclidean coordinate. If we interpret z as a time coordinate, then this equation 
describes the motion of the field in the inverted potential — Vq + ^fi'^c/)'^. Starting from an 
intial value (j){z = 0) = (pmiUai, the field will roll down the potential, pass through the 
minimum at (f> = 0, and come to rest at some field value (f) final (and then roll back towards 
the origin, but we are not interested in this motion of the field after it comes to rest at 
final)- 4>* is defined as that initial value (j)initiai for which (j) final = 0+- The solution to 
Eq.( [A.2| ) (with the boundary condition (j){z = 0) = 0=,,) is 

</,(^) = 20,^^i^ (A.3) 

where Jn is the Bessel function of the first kind of order n. The field begins with zero speed 
at 4>^,. When its speed vanishes again at the other side of the potential (at z = z^, say), 
its field value is (j)'^. It can be shown (using the recursion relations '^—^^ = Jo{z) — J2{z), 
and 2Ji(z) = z{Jq{z) + J2{z))) that this happens when J2{z = z'_^_) = 0. This gives the 
values z ~ 5.14, cj)', ~ (/)*Jo(5.14) ~ — 0.13(/>*. Since we assume and are very close 



case where the potential is fitted with an inverted potential of polynomial n (i.e., cubic for n = 3, etc). 
One shall find (and we have checked this numerically) that the thin wall approximation holds for larger and 
larger e as n is increased. 
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to each other, we get (j)^ in terms of (pj^ as 



(t>, ~ -7.7(t>+ (A.4) 
The corresponding vacuum energy difference is 

e* = \l^'^\(l>\- <t>l\- SS.SHarrier (A.5) 

where V^arrier = ^(0) — y(</>+) is the height of the barrier infront of the false vacuum 
For e > e*, the bounce becomes insensitive to the true vacuum at (p^. For e > e*, on bubble 
nucleation, the field value at the center of the bubble is 0*- It is interesting that the value 
of e* can be of the order of ten times the height of the barrier. 

The thin wall approximation has, however, broken down much earlier. To work with 
in thin wall approximation, the wall thickenss much be much less than the radius of 

the bubble {R = 35i/e), i.e. 

= SS'i/Li (A. 6) 

For the inverted quadratic potential we can calculate etw Si is given by 



5i = / dcf>y'Vi<p) - V{cf>+) 




For ~ (?!)_ ~ ibV2Vo//x, we get Si ~ vrVo///, i.e., etw = 35i;U ~ SttVo. So e can be of 
the order of Vf,a„jer (in this particular case, since the two vacua are close to V{(j)) = 0, 

Vo ~ ybarrier )■ 



B. Appendix: Evading the CPS no-go theorem 

After the initial appearance of this paper, Copeland, Padilla and Saffin produced a theorem 
in ref. [43] that analyzes the possibility of resonance tunneling in quantum field theory. 
Starting with a set of assumptions, the authors show that resonance tunneling cannot be a 
first order phase transition process in a single scalar field quantum field theory. This follows 
from the authors' argument that the "bound state" condition (which is a prerequisite for 
resonance tunneling) cannot happen in field theory. There are various ways to circumvent 
this no-go theorem (as mentioned in [43] itself). In this section we explain their proof and 
show that while the proof is valid if the all the assumptions hold, the assumptions are 
inapplicable. The very framework that the authors employ to prove the no-go theorem 
breaks down in the context of resonance tunneling. 

We first briefly explain the CPS no-go theorem. As explained in Sec. |, the occurence 
of resonance tunneling requires that the resonance condition be satisfied, i.e. 




(B.l) 
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W = f ' dX'^2m{X')[E -UeffiX')] . (B.2) 



where n = 0, 1, 2, 3, .... is related to the characteristics of the intermediate vacuum 

f A3 

'A2 

A2 and A3 being the turning points (see Sec § for derivation). Eq. (|B.l| ) is precisely the 
condition for realizing a bound state of the scalar field in the intermediate vaccum. Once 
this condition is satisfied, the transmission probability given by 

is enhanced. If the two barriers are shaped such that G = then the transmission 
probability is order unity. 

It is possible to construct an 0(4) symmetric bounce when discussing tunneling across 
a single barrier. One works entirely in the Euclidean regime (the under-the-barrier region) 
and finds the bounce solution. Analytically continuing it to the Lorentzian region gives 
the 0(3, 1) growing bubble solution. Eg. ( p. 3] ), however, describes the transition from one 
vacuum to another which lacks an entirely Euclidean description. The intermediate vacuum 
is in the Lorentzian region and contributes the oscillatory contributions of cos'^ and sin^ W. 
The Euclidean 0(4) bounce solution will not be present anymore due to these Lorentzian 
interventions from the intermediate vacuum. The usual description of the phase transition 
via 0(4) symmetric bounce is no more valid. 

Now we explain the CPS no-go theorem. Consider a 1 + 1 dimensional field theory 
for simplicity (the arguments are easily generalizable to more dimensions). The theorem 
depends on the following assumptions 

1. the bound state condition applies, i.e. W = {n + l/2)7r. 

2. it satisfies and classical equations of motion and is not the false vacuum. 

3. it obeys the boundary conditions : 0(r = 0) = (ptrue, 4>{f = 00) = 4>faise- 

4. it has zero energy. 

5. dip/dt vanishes at two different times ti and t2- 

Starting with these assumptions, the authors construct an integral I{x) given by 



I{x) 



Because of assumptions 2 and 5, dl/dx = 0. So I(at any point)= /(oo) = (since the 
field configuration at 00 is assumed to be the false vacuum, /(oo) vanishes). The energy 
functional is given by 



1 —2 




E= I d'xlV{cP) + ^l^] +^(^) ), (B.5) 
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and by assumption 4: E = 0. Combining the two results / = and E = then imphes the 
following 

That is, subject to the above assumptions, if the field configuration is sitting at the false 
vacuum at infinity, it is doomed to sit at the false vacuum every where else. 

The above no-go theorem uses several well known facts about instantons in quantum 
field theory. Because the imaginary part of the vacuum energy density (and thus the decay 
rate per unit volume) is proportional to exp(— S^;), the Euclidean action of the instanton 
must be finite. This requires the instanton to asymptote to the false vacuum. Another 
property of resonance tunneling instantons (at least in quantum mechannics) is that they 
possess a classically allowed region where the field could oscillate in real time. Because 
classical oscillation implies the existence of two turning points, CPS point out that in 
some coordinates, the field configuration describing the expanding bubble (which is also 
the instanton) should have two "spacial" slices ti and t2 over which the field is "at rest." 
(We use quotations because it is difficult to translate the analogous quantum mechanical 
problem to quantum field theory.) 

The mathematical machinery of the CPS theorem then proceeds to make the following 
argument: Since the spacial slices must extend to spacial infinity, they must extend into a 
region where the field is asymptotic to the false vacuum. Furthermore, it can be elegantly 
shown that any field configuration which obeys the equations of motion in some four- 
volume V4 is uniquely determined by its value (and derivative) on dV^. Thus since the field 
is constant on the two spacial slices ti and t2, as well constant at spacial infinity, the field 
must be simply the false vacuum solution of zero Euclidean action, and not the instanton 
we were attempting to find. Since the instanton satisfying the stated conditions does not 
exist, there is no associated tunneling event. 

The way we can evade this theorem is to point out that by continuity, the "classical 
region in which the field can oscillate" must be contained within the domain wall. This is 
simply due to the fact that the "classical" region of the potential is between the false and 
true vacua. Since the domain wall is by definition the region where the field interpolates 
between the false and true vacua, the "classical region" must be localized to a thin slice 
within the domain wall itself. This calls into question the assumption that the "spacial 
slices" ti and t2 do in fact asymptote to regions containing false vacuum. Because ti and 
t2 must be localized within the domain wall, they are by definition far (in field space) from 
the false vacuum. 

We have not shown the existence of resonance tunneling in quantum field theory beyond 
stating the necessary conditions that are to be satisfied (in Sec. ^.2]) . However we have 
shown the CPS theorem does not demonstrate the impossibility of this task. 
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